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Abstract. In this paper we shall prove that every Witten multiple zeta value of weight w > 3 attached 
to sl(4) at nonnegative integer arguments is a finite Q-linear combination of MZVs of weight w and depth 
three or less, except for the nine irregular cases where the Riemann zeta value (^{w — 2) and the double 
zeta values of weight w — 1 and depth < 3 are also needed. 

1. Introduction 

It is well-known that suitably defined zeta and i-functions and their special values often play significant 
roles in many areas of mathematics. In |21j Witten studied one variable zeta functions attached to various 
Lie algebras and related their special values to the volumes of certain moduli spaces of vector bundles of 
curves. Zagier |22| (and independently Garoufalidis) gave direct proofs that such functions at positive 
even integers are rational multiples of powers of tt. More recently Matsumoto and his collaborators 
[12), [T3l [Ml ITSl HT] defined multiple variable versions of these functions and began to investigate their 
analytical and arithmetical properties. 

Let N be the set of positive integers and No = N U {0}. For any d G N we let [d] — {1, . . . ,d) be a poset 
with the usual increasing order. By i = (ii, . . . C [d] we mean a nonempty subset of [d] as a poset. 
We say the length of i is lg(i) — £ and the weight of i is wt(i) = ii + ■ ■ ■ + ig. We define the generalized 
multiple zeta function of depth d as 

oo lg(i) 
mi,...,md = l iC[d] j = l 

For example, the Euler-Zagier multiple zeta function [U [551 US] 

C(si,...,s<i):= E mr'm^'^^---m-'- (2) 

mi>--->md>l 

corresponds to the special case that Si = unless i = {1,2, . . . ,£) for £ = 1, . . . ,d. The Mordell-Tornheim 
multiple zeta function defined by pUl) (see [HI [50]) is the case where Si = unless lg(i) = 1 or lg(i) = d. 
If we set Sj = unless i = {a,a + 1, . . . ,b) is a consecutive string of positive integers then we get exactly 
Witten multiple zeta function associated to the special hnear Lie algebra (see [T7]V 

oo 

Csi{d+i){{si,j)i<i<3<d) E n {mi+mi+i-\ hTOj)"'*'-^ (3) 

mi,...,m<j = l l<i<j<d 
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The generalized multiple zeta-functions defined by ([T]) are special cases of the functions studied by 
Essouabri [5], de Crisenoy [5], and Matsumoto ,16j. In particular we know that (d{s) has meromorphic 
continuation to the whole complex space C^**"^. However, in the form ([T]) we may have better control of its 
arithmetical properties, namely, we may be able to compute explicitly their special values at nonnegative 
integers. As usual we say a value of the Euler-Zagier multiple zeta function at positive integers a multiple 
zeta value (MZV for short) if it is finite. Our major interest is to solve the following problem. 

Main Problem. Suppose s G N^"*^^ (resp. s £ Nq and Q{s) converges. Is Cci(s) always a Q-linear 
combination of MZVs of the same weight (resp. same or lower weights) and of depth d or less? 

When d = 2 the function in llj becomes Mordell-Tornheim double zeta function. By the main result of 
[M] (see Prop.Hm know the above Main Problem has an affirmative answer for all Mordell-Tornheim 
multiple zeta functions. In this paper we will consider Cs[(4)(s) which is essentially the case when d = 3. 

Throughout the rest of the paper whenever we say some special value with positive (resp. nonnegative) 
integer arguments is expressible by MZVs we mean that the value can be expressed as a finite Q-linear 
combination of MZVs of the same (resp. same or lower) weights and the same or lower depths. Our main 
result is the following theorem which provides an affirmative answer to the above Main Problem for the 
case d — 3. 

Theorem 1.1. Suppose s S N"^ (resp. s G NqJ. // Cal^) converges then it is expressible by MZVs. 
Moreover, very Witten multiple zeta value of weight w attached to s[(4) at nonnegative integers is a Q- 
linear combination oj MZVs of weight w and depth three or less, except for the nine irregular cases defined 
by (I13|) to (I17p where the Riemann zeta value C,{w — 2) and the double zeta values of weight w — 1 and 
depth < 2 are also needed. 

The first author wishes to thank Max-Planck-Institut fiir Mathematik for providing financial support 
during his sabbatical leave when this work was done. The second author is supported by the National 
Natural Science Foundation of China, Project 10871169. 

2. Some preliminary results 

In this section we collect some useful facts which will be convenient for us to present our main result 
in later sections. 

2.1. Convergence domain. We assume that all components of s are integers and derive the necessary 
and sufficient conditions for ([T]) to converge although the conditions are still sufficient if we allow complex 
variables and take the corresponding real parts of the variables in the conditions. 
Recall that the MZV in @ converges if and only if 

E^.>^ (4) 

foi all £ — 1, . . . ,d. It is straight-forward to see the same holds for the "star" version of the multiple zeta 
function: 

C(si,...,Sd) := J2 m^''m^''---m-''. (5) 

mi >■ ■ ■ >md> 1 

Special values of ([5]) were studied in P] and [T2] . We can extend this convergence criterion easily to our 
generalized multiple zeta functions. 

Proposition 2.1. The generalized MZV 

oo lg(i) 
mi,...,md=l iC[d] j = l 
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converges if and only if for all £ — 1, . . . , d and all i — (ii, . . . , if) C [d] 

j contains at least one of 

Proof. The idea of the proof is similar to that of [24l Thm. 4]. First we observe that for any subset 
i = {ii, . . . ,ii) C [d] we have 

e 

maxjmij ■ 1 < j < £} < ruij < £raax.{mi^ ■ ^ ^ j ^ £}■ 

Hence 

oo lg(i) 

E E n(™-(^i))"'"" ^ E n(E"^^.)"''=^^-LHs, (8) 

•n-eSd m„(i)>--->m„(j)>l iC[(i] mi,...,mj = l iC[d] j=l 

where LHS means the quantity at the extreme left of the above inequalities. Observe that for each 
fixed TT G &d the power of m^(^jj in LHS is — s^(i) where i runs through all subsets of [d] whose first 
component is j. Hence the criterion ^ implies that LHS and (i(LHS) of ([5]) converges if and only if for 
each fixed vr e 

i 

E E > 

il—l i: first component is ii 

for all £ = 1, . . . ,d. Let ij = 7r^^(j) {j = 1, . . . , £). Then in the above sum i runs through all subset of [d] 
containing at least some ij {j = 1, . . . ,£). This is exactly (O, as desired. □ 

Remark 2.2. It is easy to see that similar result holds if we replace the factors '^rrii. in ^ by linear 
forms of TOi, . . . , rud with nonnegative integer coefficients. 

2.2. MZVs with arbitrary integer arguments. In this paper we will mostly consider multiple zeta 
values with nonnegative integer arguments as long as they converge. However, we will prove a more 
general result as follows since the inductive proof forces us to do so. 

Proposition 2.3. Suppose si,...,Sd £ Z. // C(si, . . . , s^) converges then it can be expressed as a Q- 
linear combination of MZVs (at positive integer arguments) of the same or lower weights and the same 
or lower depths. 

Proof. We prove the proposition by induction on the depth. When d = 1 we have nothing to prove. 
Suppose the proposition holds for all MZVs of depth d — 1. Suppose further 

si + ■ ■ ■ + Sj > j for all j — 1, ... , d, (9) 

so that C(sij • • • I Sd) converges. Assume —t = Sj < 0. Then by definition 



' m T _ 1 — 1 m T 



m*. 



c{si, ■ ■ ■ , Sd) — — 5— 1 5771 —77 ^ 

Now by the well-known formula (see, for e.g., [TTJ p. 230, Thm. 1]) 

n 

^ m* = -—(^Bt+,in+l) - Bt+m 

m— 1 

where Bt+i{x) is the Bernoulli polynomial we immediately see that ({si, . . . , Sd) is a Q-linear combination 
of MZVs of the forms Ci^i, • • • , Sj^i — s, s^+i, . . . , Sd) and ('(sii • • • i sj-i — s, Sj+i, . . . , Sd) where s = 
0, 1, . . . , —Sj + 1. All of these MZVs are easily to be shown as convergent values by (|9]) and therefore the 
proposition follows from induction assumption. □ 
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2.3. Mordell-Tornheim zeta functions. They are defined by (see [l8l [20]) 



Cmt(si 



> sd; s) 



E 



• ™d """{mi H h rud) 



(10) 



mi,...,md = l 



The main result of [24] is the following 



Proposition 2.4. ( [24) Thm. 5]) Let si, . . . , Sd and s be nonnegative integers. If at most one of them 
is equal to then the Mordell-Tornheim zeta value ^mt(si, • ■ • , s^; s) can be expressed as a (^-linear 
combination of MZVs of the same weight and depth. 

In this paper we will only need this proposition when the depth is three. 
2.4. A combinatorial lemma. The next lemma will be used heavily throughout the paper. 

Lemma 2.5. ( |24l Lemma 1]) Let r and ni, . . . ^n^ be positive integers, and let xi, .. .,Xr be non-zero 
real numbers such that a;i + • • • + a;,. 7^ 0. Then 



3 = 1 3 



3=1 \ k=l ak=0 , 
k¥=3 



k=l 
k^3 



where the multi-nomial coefficient 

(rij +Aj-l)\ 



Mi 



-1)! -TT 1 ^ . V- 

II _ and A.^^Ok. 



k=l 



k=l 



The notation means the multiple sum • • • • • • . 

ai— aj_i— Oaj + i— a^—O 



k—1 afc— 



By definition JT]) 
C3(si, • ■ • ,S7) := ^ 



3. Proof of Theorem 11.11 



1 



(11) 



^ , m^^mi'^mV' (mi + 7712)** (m2 + msY^ (mi + msY'^ (mi + m2 + ms)"'' 

ii,...,m3 — 1 -L^o\ y\ /\ /\ / 

To guarantee convergence by ([7]) we need to assume: 

i = (1) : si + S4 + se + sy > 1, 

i=(l): S2 + S4 + S5 + S7 > 1, 

i = (1) : S3 + S5 + S6 + S7 > 1, 

i = (l,2): S1 + S2 + S4 + S5 + S6 + S7 > 2, (12) 

i = (1,3) : si + S3 + S4 + S5 + S6 + S7 > 2, 

i = (2, 3) : S2 + S3 + S4 + S5 + se + S7 > 2, 

^ i = (1, 2, 3) : si + S2 + S3 + S4 + S5 + S6 + S7 > 3. 

We now use a series of reductions to prove the theorem. All of the steps will be explicitly given so that 
one may carry out the computation of pip by following them. 

Step (i). Since 2(mi +m2 +m3) = (mi +m2) + (m2 +m3) + (mi +m3) we see clearly that we can assume 



either S4 = or S5 = or sg = 0. In fact in Lemma 12.51 taking r = 3, a;i = mi + m2, X2 — m2 + m3 
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— S5 and n3 = 


: sg we get: 


(S6 + 04 + 


«5-l)! 




- I)!a4!a5! 


(S5 + 04 + 


ae-l)! 


2«7+a4+a6(s5 


- \)W-a^\ 


(S4 + 05 + 


ae-l)! 




- I)!a5!a6! 



S4 — 1 S5 — 1 

C3(S1, . . . ,S7) = ^ ^ r,,,,r°,4^J/„' "^Mr^l- , C3(S1,52,S3,S4 - 04, S5 - a5,0, Se + S7 + 04 + flg) 
a4— as— 

S4 — 1 S6 — 1 

7(3(51, S2, 53, S4 - 0^4,0, S6 - ag, S5 + S7 + 04 + ae) 

a4— ae— 

S5-I S6-1 

7(3(51, S2, 53, 0, 55 - as, S6 - ae, 54 + sy + 05 + ae). 

/r I ■ ■ -u I — II 

as— as— 

We recommend the interested reader to check the convergence of the above values by (fT2|) . The rule of 
thumb is that if we apply Lemma [2751 with each Xj a positive combination of indices then the convergence 
is automatically guaranteed. In each of the following steps we often omit this convergence checking since 
it is straight-forward in most cases. The only exception is (|26p which in fact poses the most difficulty. 

By symmetry we see that without loss of generality we only need to show that C3('Sii ■ • ■ , S5, 0, se) is 
expressible by MZVs. This is nothing but the Matsumoto's version of Witten multiple zeta function of 
depth 3 associated to the special linear Lie algebra 5l(4) (see ([3])): 

00 ^ 
C5i(4)(si, . . . , se) := ^ 



mf^mo^TO^^ (mi + m2)'**(m2 + m3)'*=(mi + m2 + 

mi,...,m3 — 1 i^o\ I \ / \ / 



Before going on we need to define the so called regular and irregular special values of C5[(4)(si) • • • > se)- 
Clearly the following special values are expressible by MZVs of mixed weights: ({0}^ means to repeat 
k times) 

/- ^^n^ un \ r ^rm h \ /-/ u n\ J C(a, & - 1) - C(a, 6), a>l,6>l, 
Csi(4)({0}3,^0,a) = CM(4)({0}4,^a) = C(a,^0) = <^ , , (13) 

[ ({a - 1) - C(a) - C(a, 1), a > 1, 6 = 1, 

C3i(4) (6, {0}4, a) = C.i(4) (0, 6, {0}3, a) = Cm(4) (0, 0, 6, 0, 0, a) = 0, b) 

^ f C(a-l,fe)-C(a,fe-l)-C(a,6), a>2,b>l, 
\ C(a-l,l)-C(«-l) + C(a)-C(a,l), a>2,5=l, 

C.i(4)({0}5,a) = C(a,0,0) = ic(a-2)-^C(«-l) + C(a), a > 3. (15) 



(14) 



By Step (ii) we will see that special values 

Csi(4)(0,0,S3,S4,0,S6) and Cs[(4)('Si, 0, 0, 0, S5, se) (16) 

are also expressible by MZVs of mixed weights. Further, if a, 6 > 2 then we have 

00 ^ 

C[(4)(0,0,0,a,6,0) = Yl 7 1 w 1 ^ (1^) 



mi ,r?i2 ,m3= 
00 / 



= E E + E + E 

m2=l \mi<m3 m\=m,-^ m\>mz/ 

<(a,6,0) +C(6,a,0) +C(a + 6,0) 



1 



(mi + m2)"-{m2 + ma)^ 



=C(^ a - 1) + C(a, 6 - 1) + C(a + 6 - 1) - C(a)C(^)- 

We call the values appearing in the nine cases from to (|17p irregular values. Otherwise 
Cs[(4)(si, . . . , sg) is called a regular value. 
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Step (ii). Taking xi = mi and X2 = m2 + in Lemma [231 we get 

SB-l 



a5=0 V "5 / 

+ ('^^'^a^^ '^^Cs[(4)(si - ai)'S2,S3,S4,0,S6 + + ai)- 



01=0 

We see clearly that we can assume either (ii.l): S5 = or (ii.2): si = and S5 > 1. Moreover, 
by taking S2 = S3 = S4 = in (fTSl) we see that Cs[(4)('Sii 0) Oi 0, S5, sg) in (fT6]l can be expressed 
by Q-linear combinations of MZVs appeared in the (jl3p and (jl4p . The argument is similar for 
Cs[(4)(0, 0, S3, S4, 0, sg) in (fT6]) . Therefore in what follows we assume that Cs[(4)('5i, • • • , sg) are 
always regular and show that they can be expressed by Q-linear combinations of MZVs of the 
same weight and depth three or less. 

Step (ii.l)- Let S5 = 0. Then we must have either si > 1 or S2 > 1 since we assume 
C0[(4)('5i) • • • ) ■54, 0, sg) = — ji — — IT? ; 7777 1 1 w 

mi,. ..,7713=1 1 z o V 

is regular. Then we may use x\ = mi and X2 = m2 in Lemma 12.51 to get 

si-l 



( S2 + Oi — l\ 

Cs[(4)(si7---)'54,0,S6) = ^ jCs[(4)('5l - ai,0,S3,S2 + S4 + ai,0,S6) 

^ /si + 02 — 1\ ^ , 

JCs[(4)(0,S2 - a2,S3,Si + S4 + 02,0,56) 



02=0 

By symmetry we only need to consider 

00 ^ 

Cs[(4)(0,S2,S3,S4,0,S6) = V sa S3 / \ wTT \ \ w (19) 

' '-^ m^mn^^ymi^ mor^ymi^ mo^ m-3,Y^ 

r?ii,...,m3=l ^ o \ 

where S2 ^ 1- But now we may take xi = mi + 7712 and X2 = m^ in Lemma 12.51 to reduce it to 
either S3 = 0, S4 > 1 or S4 = 0, S3 > 1. Then we get the following two kind of values: 

Cs[(4)(0,S2,0,S4,0,S6) = C('S6,S4,S2), Cs[(4) (0, S2, S3, 0, 0, Sg) = Cmt(s2, S4, 0; Sg), 

where S2,S3,S4,S6 > 1. But Cmt(s2, S4, 0; sg) is expressible by MZVs of the same weight and of 
depth three by Prop. [231 Case (ii.l) is proved. 

Step (ii.2). Let si = and S5 > 1. Consider 

00 

Cs[(4)(0,S2, . . . ,S6) = ^ 



mi,...,mz 



^ mf^m'^^{mi + m2Y'^{m2 + m^Y^{mi + m2 + m^^Y^ 



To guarantee convergence we must have S4 + sg > 1, S3 + S5 + sg > 1, S2 + S4 + S5 + sg > 2, 
S3 + S4 + S5 + sg > 2, and S2 + S3 + S4 + S5 + sg > 3. Taking xi = mi + m2 and X2 = m^ in 
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Lemma 12.51 we get 

^ Cs[(4)(0, S2, S3 - as, 0, S5, S4 + S6 + as) (20) 

^Cs[(4)(0,S2,0,S4 - a4,S5,S3 + Se + 04). (21) 

So we may assume that either (ii.2.1): si = 54 = 0,55, se > 1 or (ii.2.2): si = S3 = and 
S2)S4,S5 > 1, or (ii.2.3): si = S2 = S3 = and S4,S5 > 1. Here (ii.2.1) comes from ((20]) while 
(ii.2.2) and (ii.2.3) come from dH]). 



Cs[(4)(0,S2, . . . ,S6) = ^ 

03=0 
S4— 1 

+ E 

04=0 



S4 + as 
as 

ss + ^4 
a4 



Step (ii.2.1). Let si = S4 = 0, S5, sq > 1. Then we must have either S2 > 1 or S3 > 1 in 



Cs[(4)(0,S2,Ss,0,S5,S6) = ^ 



^-^ _ rrig^mg^ (m2 + ms) '^-' [mi + m2 + Jn-s)* 



mi,...,m3 



since we consider only regular values only. Thus we may put xi = 7712 and X2 = in Lemma 
I to get 



/ s + a — 1\ 

Cs[(4)(0,S2,S3,0,S5,S6) = ^ i^^ j Cs[(4) (0, S2 - 02 , 0, 0, S3 + S5 + 02 , Se) 

E*^ ^ /s2 + as ~ 1\ 
j W(4)(0,0,S3 - a3,0,S2 + S5 + a3,S6) 



a3=0 
S2-1 



yZ (^^ ^ ) C(S6, Ss + S5 + a2, S2 - 02) 

+ I S2 + as ^ S2 + S5 + as, ss - as). 



a3=0 

It is easy to see that all the triple zeta values above have the same weight. We remind the reader 
that to determine the weight of a MZV it is not enough just to add up all the components. One 
also needs to check that every component is positive. 

Step (ii.2.2). Let si = ss = and S2,S4,S5 > 1. By (fT2]) . to guarantee convergence of 

-1)^2 



Cs[(4)(0,S2,0,S4,S5,S6) = ^ 



mi,...,ni3= 



_^ (-m2)''2(mi + m2Y'^{m2 + m^Y^^{mi + m2 + ms)' 
we need to assume 



S4 + S6 > 1, S5 + S6 > 1, S4 + S5 + S6 > 2, S2 + S4 + S5 + Se > 3. 



(22) 
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Moreover, since Cs[(4)(0, S2, 0, S4, S5, sg) is regular we must have S2 + sq > 1. Putting xi = —m2, 
X2 = mi + m2 and X3 = m2 + in Lemma 12.51 we get 

Cs[(4)(0, S2,0,S4,S5,S6) (23) 
S4-1 S5-I , —1)! 
= ~1 ^T\r"l i—(-l)^'Cs[(4) (0,0,0,54 - 04,55 - 05, S2 + S6 + 04 + 05) 

^-^ So — 1 !a4!a5! ' 

04=005=0 ^ ^ / 4 

+ —, — ^rrrn — (-i)"'Cs[(4)(o,s2 - 02,0,0,55 - 05,54 + se + 02 + 05) 

02=005=0 

S2-I S4-I / ^ _|_ + 04 — 1)! 

+ yi — ^T\rn — j— (-i)"'Cs[(4)(o,s2 - 02,0,54 - 04,0,55 + S6 + 02 + 04)- 

02=0 04=0 

We point out that since we have used x\ = —m2 we need to check the convergence of ah the 
above three kinds of values under the assumption ()22p even though the checking itself is trivial. 

Returning to the reduction of ()23p we see that the last two sums are expressible by MZVs so 
we only need to consider those values appearing in the first sum, namely, those of the form in 
the next case. 

Step (ii.2.3). Let si = S2 = S3 = and 54,55 > 1. Then we must have se > 1 since we only 
consider regular values. To guarantee convergence of 

00 ^ 

Cs[(4)(0,0,0,54,55,56) = V 7 ■ —7 ■ 7— ■ ■ — (24) 

mi,...,ni3=l ^ ' ^ ' ^ ' 

we need to assume 54 + 56 > 1, 55 + 53 > 1, and 54 + 55 + 55 > 3. If either 54 > 1 or 55 > 1 
then we may further assume that 54 > 1 by change of index mi ^ 7773 in (|24p. We may set 
x\ = —{m2 + 7773) and X2 = mi + 7772 + 7/73 in Lemma 12.51 to get 

-2 



f s + a l\ 

Cs[(4)(0,0,0,S4,S5,56) = ^ 1^ ^ J (-l)''*Cs[(4)(s6 + 05,0,0,54,55 - 05,0) (25) 

+ ^( ^ ^ I (-l)''^Cs[(4)(s5 +06,0,0,54,0,56 - ae) 

06=0 V «6 / 



^^56 + 55 2\ + 1^0,0,54,0,1) 

-Cl;^,\(56 + S5- 1,0, 0,54, 1,0)), 



■'s[(4)^ 

where C}-^^ is the partial sum of ([T]) when each index 777^^ goes from 1 to N . Observe that all 
the values in the first two lines above are convergent and the second line is already expressible 
by MZVs of the same weight. So we are reduced to consider the following two kinds of values: 

(A) . Cs[(4) (■si, 0, 0, 54, t, 0) (with 5i,54 > 1 and t > 1) from first line above and 

(B) . The limit ([26]) (with 5 := 55 + 56 - 1 > 1). 

We now need to divide into two subcases to compute ([25|) (leading to (A)) and ([26|) (leading to 
(B)): (ii.2.3.1) 5i = 52 = 53 = 0, either 54 > 1 or 55 > 1, and (ii.2.3. 2) 5i = 52 = 53 = and 
54 = 55 = 1. Before treating the two cases separately we first deform the expression inside the 
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limit (I26p as follows: 

C^J2)(s> 0, 0, S4, 0, 1) - Cij2)(«' 0, 0, S4, 1, 0) 

N /mi+rn2+N mi+m2\ ^ N /m2+N m2 \ ^ 

= y y - y - y y _y I 

mi,m,2=l \ n=l n=l / mi,m2=l \ n=l n=l/ ^ 

A'' mi+m2 N mi+m2+N ^ 

(mi + m2)*-*ri ~'~ mf (mi + m2)'*'in ^ 

mi,m2=l n=l+m2 ^ rrii ,m,2=l n=l+m2+A'' ^ 

Step (ii.2.3.1). Let si = S2 = S3 = 0, sg > 1, either S4 > 1 or 55 > 1. Without loss of generality 



assume S4 > 1 as we have remarked after ([2 
(A). Since si > l,S4,t > 1 we have 

00 

Cs[(4)('Sl,0,0,S4,t,0) = ^ 



_^ m^^ (mi + m2)*'' (m2 + ma)* 

y y I 

Z__/ Z__/ m,^^(m,i -I- mo"!**? 



mi,m2,rri3 
00 



mf^ (mi + m2)**?^' 

mi,m2=l ri>»Ti2 ^ ' 



00 ^ 

--C{t)C{s4, Si) - Cmt{si, t;s4) - y yZ [ \in 

mi,m2=ln<m2 ^ ^ ' 
00 

<(0C(s4,Sl) - CMT(si,i;S4) - X] 



mi,n,m3 



^ m^^ (mi + n + 7x1-3)^^11^ 



by setting m2 = n + ma. But the last sum is just Cmt(si, t, 0; S4) which is expressible by MZVs 
of the same weight and of depth three by Prop. 12. 4[ 

(B). The limit (|26p . This has been reduced to (|27p with s > 1 and 54 > 2. Then the second 
sum in ([27|) is bounded by 

1 ^ 1 
— - y = 0(log2 iV/iV) ^ 0, asiV^oo. 

nil ,m2=l 

By setting m2 = ma + n in the first sum in (j27p we get 
^lirn^ (Cij!;)(s, 0, 0, S4, 0, 1) - d^^lis, 0, 0, S4, 1, 0)) 

CO OD / m\+m2 "12 \ 1 

= _yy y _y i 

mi n=l \ n=l n=l/ i \ ^' 

CO /mi — 1 m2 \ ^ 

~ + X] "^X^ X^j mf (mi + m2)'*4n ^'^^^ 

mi,m2=l \ n=l n=mi m\<n<m\+m2 n=mi+m2 n=l/ ^ ' 

00 m2 — 1 _o _i 

=Cmt(s, 1; S4) - (C(S4, S, 1) + C(S4, S + 1) + C(S4, 1, S) + C(S4 + 1, «)) + Yl \ ^ )S4 ' 

mi,m2=l n=l ^ ^ 

In the last sum setting m2 = n + ma then we get 

CO m2 — 1 1 °° 1 

Y] ~rf [ = ~T( \ 1 \s. = CMT(g,i,0;s4) 

^-^ ^-^ mi(mi+m2) *n ^-^ m?(mi + ma + n)**n 

mi,m2=l n=l ^ n,mi,m3=l ^ ' 
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which is expressible by MZVs of the same weight and of depth three by Prop. 12.41 

Step (ii.2.3.2). Let si = S2 = = and 54 = S5 = 1. Then we must have se > 2 in order to 
have convergent values. Moreover, ()25p is vacuous so we don't need to consider (A). 
(B). Observe that the second term in ()27|) is bounded by 

1 ^ 1 

— - ^^ = Oi\og\N)/N)^0, as TV ^00 

since s > 2. So 

TV mi+r?i2 -. 

lim Mj^Us, 0,0, 1,0,1) -CS)(^, 0,0, 1,1,0)) =- lim ^ — 

mi ,m2=l n=l+m2 

is expressible by MVZs by Lemma 14.21 (with t = 1) in ^ 

This finishes the proof of the theorem. □ 

4. Regularized MZVs and a technical lemma 

We need to derive several formulas concerning MVZs. They are best understood as conse- 
quences of the stuffle (also called quasi-shuffle or harmonic product) relations for the MZVs 
including the divergent ones suitably regularized. 

By [To] one may define a regularized value C(s) (denoted by Z™(T) in [10]) if s = (si, . . . , s^) 
with si = 1 such that the following stuffle relation holds: 

C(si)C(s2) = C(si * S2) (29) 
where the the product * is the stuffle product (see [9] or [ini §1]). 
Lemma 4.1. Let s and t be two positive integers greater than 1. Then 



S t 

as)m = Y.(^~^^'')at+a,s-a)+Y,{'~\'^\{s+b,t-b). 

n—n \ ^ J h—n ^ ^ 



(30) 



a=0 ^ ' 6=0 

Further 

C(s + 1) = ^C(l + a,s-«)- (31) 

a=l 

Proof. Equation ([30]) is the famous Euler's decomposition formula [7] (or see [3]). To derive (j3T]) 
we use the regularized MZVs as follows. For any positive integer A'^ we have by Lemma 12.51 

N N-m ^ N N-m / ^ s-1 ^ \ 

^ ^ mn" ^ ^ ( rn(m + nY ~^ ^ n^-^-im + nY+'^ ] ' 

m=l n=l m=l n=l \ ^ ^ a=0 ^ ' / 

Taking the regularization we get 

C(l)C(s) = C(l, s) + 1) + 5] C(l + a, 5 - o). 

a=l 

On the other hand by stuffle (f29]) 

C(i)C(s) = C(i,s) + C(s,i) + C(5 + i). 

Equation ([3T]) follows immediately. □ 
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Lemma 4.2. For all positive integers s and t such that s + 1 > 3 the value 

oo mi+m2 ^ 

y y — - — w (32) 

mi ,m2=l n=l+m2 

is expressible by MVZs of weight s + t + 1 and depth three or less. 
Proof. Let's consider the partial sum 

m-j / ms m3— mi\ 



= E E- E 



l<mi<m3<N n=l+m3—mi ^ l<mi<m'j,<N \n=l n=l / i o 



whose hmit when N ^ oo \s clearly (j32p . Using regularization and playing similar trick of 
breaking the sum as we did in ([28]) we find that 



-M E E-Eh;^. 

I l<mi<mz<N \n=l n=l / ^ ^ ) 

<(l,t,s) + C(l,s + t) + C(l,5,i) + C(i + l,s)-Reg<^ V V ^^\. (33) 

I z — z — ^ m^m^ri I 



ms—mi 



K l<mi<m-j,<N n=l 

Now the untreated sum can be deformed as follows: taking m2 = — nii we get 

"12 — 1 
l<mi<mi+m2<N n 



"12-1 



=1 1*^ ' l<mi<mi+"i2<Ar 1 2V 1 I 



= E E E . ' , +CS(^,t;l) (34) 

"11=1 n=l "3=1 

by setting 713 = 777-2 — n. Taking xi = nii and X2 = tt- in Lemma 12.51 we see that ()34p becomes 

AT A-mi-1 Af-"ii-" C s-1 , 

E E E ^ 



mi=l "=1 "3=1 ^. a=0 

t-1 



a J ml "(7771 + 77)*+"(?77i + 77,3 + n) 



^ ^ ^ 6 y 77* ^(7771 + 77)*+^(?77i +773 + 77) 



Taking regularization and combining with ([33]) yields 



m =C(1, t, s) + C(l, s + t) + C(l, s, t) + C(t + l,s)- Cmt(s, t; 1) 



I. a=0 ^ ' 6=0 

<(i,s + t) + C(t + i,s)-CMT(s,t;i) 

r s~i /, , -,\ t-i 



EC^r')^'(i'*+"'^-«)+^C^6"')^'(i'^+^'*-'4- 

a=l ^ ^ 6=1 ^ ^ I 
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Applying the stuffle relations (|29l) we get 

^ =C(l)C(s + t)-C{s + t + l)- C{s + t,l)+ ({t + 1, s) 




1 



t-i 



s + b 
b 



CMT(s,i; 1) 



C{s + b,t-b) 



a=l 



t-1 
b=l 



s + b-1 
b 



C{t + a,s-a) + Y, 

' 6=1 
|C(t + a + 1, s - a) + C(t + a, 1, s - a) 

+ C(t + a, s - a + 1) + C(t + a, s - a, 1)| 
{C(s + 6 + 1, t - 6) + C(s + ^ 1, t - 6) 

+ C(s + 6, t - b + 1) + C(s + 6, t - 6, 1)} . 

We can cancel all the terms involving (^(1) in the two cases t > 1 and t = 1 by applying (|30p and 
pip of Lemma 14.11 respectively. This finishes the proof of the lemma. □ 

5. Some examples 

Using Maple we have verified all the formulas on [17, p. 1502] by our general approach. We 
can compute any special values of C3('Si, ...,57) and in particular Cs[(4)('5i5 • • • , se) at nonnegative 
integers whenever they satisfy the convergence conditions given by (jl2p . For example, by Maple 
computation there are 32 weight four convergent Cb[(4)('515 • • • > ■sg) values at nonnegative integers, 
with 15 distinct ones as follows: 



c 
c 
c 
c 
c 
c 
c 

c 
c 

Ca 

c 



iC(2) - iC(3) + |C(2)^ Ca((4)(0, 0, 0, 2, 2, 0) = 3C(3) - 



)(0,0,0,0,0,4 

) (0, 0, 0, 0, 1, 3) = Ca((4) (0, 0, 0, 1, 0, 3) = C(2) - C(3) - U^^? ^ 

) (0,0, 0,0, 2, 2) =Ca[(4)(0,0,0,2,0,2) = ((3) - ^C(2)', 

)(1,0,0,0,1,2) =Ca[(4)(0,0,l,l,0,2) =C(3)-iC(2)^ 

) (1, 0, 0, 0, 2, 1) = Ca((4) (0, 0, 1, 2, 0, 1) = 2C(3) - \C,(2f , 

) (0,0, 1,0,0,3) =Ca[(4)(0,l,0,0,0,3) =Ca[(4)(l,0,0,0,0,3) = 2C(3) - C(2) - ^C(2)^ 

)(o,i,o,i,i,i) = ^C(2)^ Ca[(4)(i,i,i,o,o,i) = fc(2)^ 

)(0,0,0,1,1,2) = Ji5C(2)^ Ca[(4)(l,0,l,l,l,0) = iIC(2)^ 
)(1,0,1,0,0,2) =Ca((4)(0, 1,1,0,0,2) =Ca[(4)(l,l,0,0,0,2) =|C(2)^ 

)(1,0,0,1,1,1) =Ca((4)(l, 0,0, 1,2,0) =Ca[(4)(0, 0,1, 2,1,0) =Cal(4)(0, 0,1, 1,1,1) = \C(2f , 
)(0,1,0,1,0,2) =Ca[(4)(l,0,0,l,0,2) =Ca[(4)(0,l,0,0,l,2) =Cat(4)(0,0,l,0,l,2) = |C(2)2, 
)(1,1,0,0,1,1) =Ca[(4)(l, 0,1, 1,0,1) =Ca[(4)(l, 0,1, 0,1,1) =Ca[(4)(0, 1,1, 1,0,1) = §((2)^. 

Notice that in the first seven values the weights are mixed which correspond to the nine cases 
of irregular values (113p to (I17p . In higher weight cases if we only consider regular (pure weight) 
values then we have the following examples. 

C,((4)(l, 1,0, 1, 1, 1) = C[(4)(0, 1, 1, 1, 1, 1) = |C(5) - C(2)C(3) = .6150150376 . . . , 
C,[(4)(l, 0, 1, 1, 1, 1) = -|C(5) + C(2)C(3) = .4219127176 . . . , 
Q(4)(l, 1, 1, 1, 1, 1) = -§^Q{^f + 2C(3)2 = .2617453537 . . . , 

C3(l, 1, 1, 1, 1, 1, 1) = iCa[(4)(l, 1, 1, 1, 1, 2) = f C(7) - iC(2)C(5) = .08840016918 . . . , 
Cs[(4)(2,2,2,2,2,2) = 8^C(2)' = .0083233212. 
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The last value of the above agrees with [12^ (4.29)] and the first value in [U Table 4] by the 
relation 

Cs[(4)(2m,2m,2m,2m,2m,2m) = 12-2"'Csi{4)(2m), 

where Cs((4) the LHS is the multiple variable function being studied in the current paper and 
Cs((4) on the RHS is the original single variable Witten zeta function (see [171 Prop. 2.1]). As 
two intriguing computation we have 

C,((4) (1, 2, 3, 3, 2, 1) = C(2)C(8, 2) + iiiC(2)^ - iC(2)C(5)2 - f C(3)C(9) - 35C(5)C(7) 

-2C(7)C(2)C(3) + f C(10, 2) = .0129650292 . . . , 

Cm(4) (3, 2,1,1, 2, 3) = 10C(2)C(8, 2) - ^C{2f - 6C(2)C(5)2 + 44C(3)C(9) + 40C(5)C(7) 

-20C(7)C(2)C(3) - 22C(10, 2) = .0056078053 .... 

We have also verified all the above computation numerically with Maple and EZface [3]. 

We conclude our paper with a conjecture. Let M.ZV{w, < I) be the Q- vector space generated 
by MZVs of weight w and depth < /. Then by Theorem 11.11 the space generated by special 
values of Cs[(4) of weight w > 3 over Q is included in the space A4ZV{w, < 3) + J^ZVi^w — 1, < 
2) + A4ZV{w — 2, 1) (conjecturally the sum is direct). 

Conjecture 5.1. The space generated by special values o/Cs[(4) of weight w > 3 at nonnegative 
integers over Q is 

MZV{w, < 3) e MZV{w - 1, < 2) e MZV{w - 2, 1). 
Using the MZV table [2] we have verified the conjecture for all weights up to 12. 
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